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We analyze quantum state estimation for finite samples based on symmetry information. The 
used measurement concept compares an unknown qubit to a reference state. We describe explicitly 
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I. INTRODUCTION 

The concept of a measurement Hes at the heart of quan- 
tum mechanics. For a long time this has been a topic of 
fundamental discussion p]. Quantum information sci- 
ence [2 has transformed it to a basic technical issue, at 
the core of many applications |3]. Of particular impor- 
tance are quantum state measurements, since they can 
reveal complete information about a quantum system, or 
more precisely about its preparation process jH H]. In 
this sense such measurement schemes attach a certain 
operational meaning to a quantum state. 

Full tomography j6] of an unknown state requires in- 
finitely many identical copies of a quantum system. The 
no-cloning theorem j7], however, prevents to produce 
them from a single system. Hence we always encounter 
the typical situation of an estimation approach |8] |9]: It 
is only possible to collect a limited amount of informa- 
tion from a finite sample of identical systems. Using the 
corresponding finite set of measurement data, we then 
have to deliver an estimate for the underlying quantum 
state. The benchmark for the quality of such a proce- 
dure is the optimal average fidelity of a joint measure- 
ment jlOl [TT] [T2] on all systems of the finite sample. On 
the other hand, the needed physical principles for such 
collective measurements are not simple to realize. 

Hence this optimal approach has been approximated 
by various local measurement schemes |13l [T4| [15] [T6| llTl 
[TS] combined with suitable classical communication. The 
physics of the related measurements on single systems has 
a rather simple operational meaning, while the obtained 
estimation fidelities can be close to the optimal ones. 

The present contribution aims at the estimation of a 
pure qubit of which only a limited number N of iden- 
tically prepared copies is available. Our aim is not to 
optimize a certain fidelity. We rather examine how a 
basic measurement concept can be used for estimation of 
quantum states. We investigate measurements that com- 
pare the unknown copies, one by one, to a known refer- 
ence state, a "quantum ruler". For these comparisons we 
restrict ourselves to measurements that yield symmetry 
information only. Such measurements have been used in 
|19) for state comparison, and based on this work applied 
to programmable unambiguous state discrimination |2i)j. 
Moreover these symmetry measurements are a simple re- 
alization of a quantum "multimeter" (21], where the ruler 



state is "quantum software" that defines what measure- 
ment is performed. 

II. BASIC CONCEPTS 
A. Symmetry measurements 

To be more specific, our comparison is done by sym- 
metry measurements. In contrast to a spin measurement 
such a comparison is close to classical concepts of mea- 
surement: We determine the agreement between the un- 
known system and a known ruler. Therefore, it will be 
interesting to see how successful this ansatz can be in the 
quantum domain. 

In order to define this in more detail, we write the 
unknown qubit state 

|^)=cos^|0) + e'^sin^|l) (1) 

in a computational basis {|0) , |1)}. It is assumed that N 
identically prepared copies of this qubit are available, i. e. 
they all stem from the same but unknown preparation 
process. Each copy is compared to a reference state 

M = cos^|0) + e^'^"sin^|l), (2) 

with 1/ = 1,...,7V. Hence for each measurement step v 
the corresponding two-particle system is in the product 
state lip) |r^). Our measurement concept now consists in 
determining the symmetry of this state with respect to 
particle exchange. The results will be dichotomic: either 
we obtain a "symmetric (s)" or an "antisymmetric (a)" 
signature. 

The probability of finding the product state in the anti- 
symmetric subspace, spanned by the antisymmetric state 
l*-) = 75(|01)-|10)),reads 

PaiW , |r.)) = 

= - (1 — cos cos t?^ — cos (0 — (jCi,) sin0sini?,y) . (3) 

We denote the complementary probability, i. e. the prob- 
ability of finding the product state in the symmetric sub- 
space, by Ps(IV') , K)) = 1 -Pa{\i^) , \ri.)). 
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Before discussing how to utilize this kind of measure- 
ment for quantum state estimation, we mention a simple 
realization by linear optics [22] ■ Let us assume, that the 
unknown and the reference qubit are prepared as super- 
positions of single photon polarizations. These photons 
impinge on a 50:50 beam splitter that is perfectly bal- 
anced concerning the two polarizations in each of its in- 
put ports. That is, the unknown photon enters one input 
port, the reference photon the other one. The probabil- 
ity of finding a single photon with arbitrary polarization 
in each output port is then equivalent to Eq. ([3|. We 
therefore can measure it with efficient detectors that sim- 
ply distinguish the vacuum from any finite number of 
photons. This measurement scheme has already been ex- 
perimentally utilized for measuring the overlap of quan- 
tum states |23] . 



B. Estimation method 

We perform symmetry measurements on N copies of 
the qubit \ip) and the corresponding reference states. For 
each single reference \r,j) we obtain a measurement result 
a^, which is either a, for the antisymmetric subspace, or 
s, for the symmetric subspace, respectively. According 
to Bayes' rule |241 [25] we can calculate the probability of 
this sequence S^r = {a^, \r^) ;v & [1, N]} of results and 
references, and use it as a likelihood function jS] |26] 

JV 

Ln^IIpc,AW,K)) (4) 
I/=l 

for the unknown qubit Note that Ljsi is given as a 
function of the parameters 9 and </>, according to Eq. ([ij . 
Their values at the maximum oi Ln determine the esti- 
mated state IV"^*) for one specific measurement sequence 

Sjv EH. 

So far we still have no rule how to choose the reference 
states. The first reference |ri) is always arbitrary, since 
we have no information about the state \^fj) before the 
first measurement. We therefore always choose |ri) = |0). 
But at a later measurement stage we already have the 
results of the preceding measurements. Then an adaptive 
strategy for choosing the new references using the already 
obtained data is possible and will be described in the next 
section. After fixing the next reference according to a 
chosen adaption strategy, we perform the next symmetry 
measurement. Then we repeat this cycle until all copies 
of the unknown state are used up. 

To quantify how well the unknown state \ip) and its 
final estimate j-iAjv*) coincide, we apply the corresponding 
fidelity Fat = |('(/'^*|7/')r. This value of course depends on 
|-0) as well as on the probabilistic outcomes {a,^} in the 
measurement sequence Sat. To arrive at a measure for 
the quality of the used adaptive measurement strategy 
we have to average over all unknown states \ip) on the 
Bloch sphere, and all measurement outcomes. This leads 
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Figure 1: Fidelities from symmetry measurements. The 
achieved mean estimation fideUty (Fm), Eq. ([sj, is plotted 
over the number A'' of performed measurements, i. e. the num- 
ber of copies of the unknown qubit. Using an adaptive method 
(symbols ▼) gives fidelities which are up to 5% above those ob- 
tained with random references (symbols ■). Note that these 
numerically simulated data points have a small error due to 
the finite amount of Monte Carlo runs. However, this error is 
smaller than the size of the symbols, as can be seen for the 
A'^ = 1 data points, which should in principle coincide 

to the mean expected estimation fidelity 

{Fn) = I {\{rN'\i^)\')^jn^, (5) 

where the integration J dfl^ = J^^ d4> Jq d9 sin de- 
notes a Bloch sphere average and is the average 
over all measurement sequences possible within the cho- 
sen adaption strategy. 

III. ADAPTION STRATEGY 

The simplest way to construct an estimation scheme, 
which achieves equal average fidelities for all possible un- 
known states, is to use equally distributed random refer- 
ence states. Of course this is not an adaption, but a lower 
bound to which adaptive schemes can be compared in or- 
der to assess their quality. To calculate the corresponding 
fidelity, Eq. we use Monte Carlo techniques. That 
means the result of a measurement is simulated by a ran- 
dom number generator, in such a way that the results 
follow the distribution Eq. ([3|. From these randomized 
results we calculate F/v, and by averaging over 10^ runs 
we arrive at a reliable value for (Fn), Eq. ([s}. The results 
obtained with random references are depicted in Fig.jl] 

As one can see, the average fidelity (F/v) grows with 
the length N of the measurement sequence. Hence in 
principle an estimation based on a symmetry compari- 
son with reference states works. However, the average 
fidelities are quite low. Therefore, we ask whether this 
quantum estimation method allows for adaption, which 
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Figure 2: Estimated states on the Block sphere. On the left 
hand side we depict the Bloch sphere distribution of estimated 
states for 10'' different runs and three sizes A'^ of the finite 
ensemble. Each point {On^ , (j}'§'^) corresponds to one simu- 
lated measurement run, where the prepared, unknown state 
has been chosen randomly. In order to see the convergence, 
the coordinate frame has been rotated after the estimation 
such that the unknown state would lie on the north pole, i. e. 
cosO = 1. On the right hand side the same data sets are 
shown as histograms for cos 6. The convergence of the esti- 
mated states towards the correct state can be clearly seen. 

means choosing appropriate references \r^) from one mea- 
surement step to the next. 

Our goal is to maximize the fidehty of the finally esti- 
mated state. We can directly use this as a quality mea- 
sure for finding adapted reference states. An optimal 
way to do this would be to calculate the fidelity of the 
last estimated state and maximize it with respect to the 
reference state we want to adapt. In order to keep the 
computational effort reasonable, we restrict ourselves to 
the maximization of the fidelity F^+i after the next mea- 
surement. 

To explain this in detail, let us assume that we have al- 
ready performed v measurements. Then we can calculate 
an expected fidelity for the next measurement step ly + 1 
before we perform it. Here "expected" means to average 
over the possible measurement outcomes. Such a quan- 
tity is a function of the next reference state \r,y+i). To 
get the best possible fidelity after the next measurement 
we have to choose a reference that maximizes it. 

To quantify this expected fidelity, we have to determine 
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Figure 3: Reference states on the Bloch sphere. Here we show 
distributions of reference states, Eq. (|2|, over the Bloch sphere 
and corresponding histograms, again rotated as described in 
the caption of Fig. [2] Due to this rotation the first reference 
states |ri) = |0) (in the calculational basis) are randomly dis- 
tributed (a) with respect to the unknown state. It can be 
seen, that for higher numbers of measurements the references 
tend to be close to the prepared state, but they do not become 
completely parallel. Also, in comparison to Fig.[2]the conver- 
gence is weaker. The azimuthal asymmetry in (b) is due to 
the small number of possible adaptions for the reference state 
jrs). After four steps we can obtain one of 2* possible sets of 
measurement data. Therefore, there are only 16 possibilities 
for the orientation of the reference Irs) relative to |ri). As can 
be seen, this breaks the equipartition over the 1^95 coordinate. 

the states I'lp'^^^i a) and \ipu%i s) that we would estimate 
if we got the outcomes a or s, respectively. Note that 
these states themselves depend on |ri,+i). Then we av- 
erage over the corresponding fidelities alV*)! a-nd 

KV'S^+i sIV")!^- Moreover, since we do not know ji/;), we 
have to average over it too, and hence arrive at the con- 
dition 

+p, (|r,+i) , 1^)) • |(C^i,.k)r } — max. (6) 

for the next reference state |ri,_|_i). 

Note that we have neglected the information on the dis- 
tribution of 1-0), we had already obtained in the preceding 
steps, namely the likelihood Li^, Eq. ([i}. In other words, 
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the V measurement outcomes and their corresponding 
reference states, i. e. the sequence Si,, enter impHcitly via 
the adaption of the reference states but not expHcitly in 
a weight factor for dfi^. We gain the advantage that all 
expressions to be integrated are known functions of the 
involved quantum states. Hence we can easily perform 
the integration and arrive at an analytical expression de- 
pending on l^jy+i). Nevertheless, here might be room for 
further improvements of a quantum adaption. To assess 
this adaption strategy, we again simulate the correspond- 
ing estimation scheme, and compare it to the random ref- 
erence strategy. As one can see from Fig.jl] application 
of the adaptive scheme gives significantly better fidelities 
than a random selection of reference states. 

To further elucidate the convergence of the adaptive 
scheme, the evolution of estimated states on the Bloch 
sphere is depicted in Fig.|2] As expected for an adap- 
tion, the estimated states accumulate near the unknown 
state. The simulation shows, that in doing so, they fol- 
low, within the numerical errors, an exponential distri- 
bution. 

Another interesting question is the distribution of ref- 
erence states. If they would prefer a certain alignment 
relative to the unknown state, e. g. parallel to them, an 
experimentalist would not have to use the quite involved 
adaption procedure described above, but could simply 
aim for this alignment. As can be seen from Fig.|3] this 
is not the case. The reference states have a clear tendency 
to be on the same hemisphere as the unknown states, but 
they show a rather broad distribution and do not become 
completely parallel. 

IV. CONCLUSIONS 

Symmetry measurements have been discussed in the 
context of quantum state estimation. This measurement 
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